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Abstract 



Net-baryon, net-charge and net-strangeness number fluctuations in high energy heavy-ion collisions are discussed within the frame- 
5 work of a hadron resonance gas (HRG) model. Ratios of the conserved number susceptibilities calculated in HRG are being com- 
^ pared to the corresponding experimental measurements to extract information about the freeze-out condition and the phase structure 



< 



of systems with strong interactions. We emphasize the importance of considering the actual experimental acceptances in terms of 
kinematics (pseudorapidity (77) and transverse momentum {pr)), the detected charge state and effect of resonance decays before 
comparisons are made to the theoretical calculations. In this work, based on HRG model, we report that the net-baryon number 
^0 fluctuations are least affected by experimental acceptances compared to the net-charge and net-strangeness number fluctuations. 
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1. Introduction 

Measurement of the moments of distribution for conserved 
quantities like net-baryon, net-charge, and net-strangeness 
number for systems undergoing strong interactions as in 
high energy heavy-ion collisions, have recently provided rich 
physics insights Q S B EL 0, S H H, E3] ■ The most cru- 
cial realization is that, the product of moments of the conserved 
number distributions are measurable experimentally and can be 
linked to susceptibilities (x) computed in Quantum Choromo- 
dynamic (QCD) based calculations [lLI^U- For example, Scr 
- X^ Ix^ an d K(jl - X Ix^i where <x is the standard de- 
viation, S is the skewness, k is the kurtosis of the measured 
conserved number distribution, x^ are the n th order theoreti- 
cally calculated susceptibilities associated with these conserved 
numbers. Such a connection between theory and high energy 
heavy-ion collision experiment has led to furthering our under- 
standing about the freeze-out conditions J2L |4|], details of the 
quark-hadron transition [l] |8J and plays a crucial role for the 
search of possible QCD critical point in the QCD phase dia- 
gram [51]. In all such physics cases there is a need to establish a 
reference point for the measurements. Computing these quan- 
tities within the framework of a hadron resonance gas (HRG) 
model lllll provides such a reference for both experimental data 
and QCD based calculations. 

The experimental measurements have limitations, they are 
usually available for a fraction of the total kinematic phase 
space, due to the finite detector geometries and can detect 
only certain species of the produced particles. For example, 
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measurements related to net-baryon number distribution is re- 
stricted by the kinematic range in pj where their identification 
is possible. In addition baryons like neutron are not commonly 
measured in most of the high energy heavy-ion collision exper- 
iments. While for the net-charge number distribution, the limi- 
tation is usually in kinematic range available in r\ and the details 
of how contribution from different charge states and resonances 
are dealt with in the measurements. The kinematic acceptance 
in a typical high energy heavy-ion collision experiment like 
STAR (H and PHENIX fH at the Relativistic Heavy-Ion Col- 
lider facility (RHIC) for net-charge multiplicity distributions 
are: |j/| < 0.5, 0.2 < p T < 2.0 GeV/c and \tj\ < 0.35, 0.3 < p T < 
1.0 GeV/c, respectively. While for net-baryon number and net- 
strangeness related studies carried out in the STAR experiment, 
within I77I < 0.5, is through the measurement of net-protons and 
net-kaons in the range of 0.4 < pj < 0.8 GeV/c and 0.2 < pj < 
2.0 GeV/c, respectively JKH. 

The main goal of this paper is to demonstrate using the HRG 
model (discussed in next section), the effect of the above ex- 
perimental limitations on the physics observables^V^ 2 ' and 
X^ IX (T> ■ Our model based study clearly shows that the value 
of the observables related to net-charge and net-strangeness 
strongly depends on the experimental kinematic and charge 
state acceptances. Where charge state could be electric charge 
02) = 1 or higher for net-charge measurements and strangeness 
number (S) = 1 or higher for net-strangeness measurements. In 
contrast, the net-baryon number studies are found to be mini- 
mally affected by these experimental limitations. 

The paper is organized as follows. In the Section [2] we will 
discuss the HRG model used in this study. In the Sectionf3j the 
results for the observable x^/X^ and ;^ (4 V;f (2) are presented 
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for different kinematic acceptances, charge states, and reso- 
nance decay contributions. We also provide a table listing the 
values of these observable for typical experimental conditions 
as encountered in STAR and PHENIX experiments at RHIC 
and ALICE experiment at the Large Hadron Collider (LHC) 
Facility. Finally in Section [4] we summarize our findings and 
mention about the implications of this work to the current ex- 
perimental measurements in high energy heavy-ion collisions. 

2. Hadron Resonance Gas Model 

In the HRG model, we include all the relevant degrees of 
freedom of the confined, strongly interacting matter and also 
implicitly takes into account the interactions that result in reso- 
nance formation yfl. The logarithmic of the partition function 
(Z) in the HRG model within the ambit of a grand canonical 
ensemble framework is given as, 



\nZ(T,p B ,p Q ,p s ) = 2_jlnZi(T,fi B ,fiQ,Hs) 

B 

+ ^lnZi(T,jj Q ,fi s ) , 



(1) 



M 

where, lnZi(T, V,p) = ±0 J<9 3 /cln{l + exp(p - E)/t], T is 
the temperature, V is the volume of the system, p t is the chem- 
ical potential and gj is the degeneracy factor of the i' h particle. 
The total chemical potential p = Bp B + Q^q + Sps, where B, Q 
and S are the baryon, electric charge and strangeness number 
of the i' h particle, with corresponding chemical potentials p B , 
jj.q andyUj, respectively. The '+' and '-' signs are forbaryons 
and mesons respectively. 

The thermodynamic pressure (P) can then be deduced for the 
limit of large volume as 



P 

T A VT 



*- InZi = ±^3 Jd 3 kln (l ± expQi - E)/t}.(2) 



The n th order generalized susceptibility for baryons can be 
expressed as, 
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and for mesons, 
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The factor X represents either B, Q or S of the i' h particle de- 
pending on whether the computed Xx represents net-baryon, 
net-electric charge or net-strangeness susceptibility. Here 
d 3 k = pr -\Jpt 2 + tw 2 cosh 77 dpj dr\ d<p and E - mj cosh 77 = 
■\Jpt 2 + tn 2 cos h 77. Where m is the mass of the hadron and cf> is 
the azimuthal angle of the produced hadron. The experimental 



acceptances can be incorporated by considering the appropriate 
integration ranges in 77, pj, <p and charge states by considering 
the values of \X\. The total generalized susceptibilities will then 

be the sum of the contribution from baryons and mesons as, 

(«) («) , («) 

V = V -4- V 

X* Xx,baryon T Ax,meson- 

In order to make the connection with the experiments, the 
beam energy dependence of p and T parameters of the HRG 
model needs to be provided. These are extracted from a sta- 
tistical thermal model description of measured particle yields 
in the experiment at various y^NN Il4l Il5l llol . This is fol- 
lowed by the parameterization of p B and T as a function of 
V^nn Jl4 Il6ll . The ps dependence of the temperature is given 
as T(p B ) = a - bp\ - cp\ with a = 0.166 + 0.002 GeV, 
b = 0.139 + 0.016 GeV 1 , and c = 0.053 ± 0.021 GeV 3 . The 
V^nn dependence of p B is given as p B ( ^$nn) - 1+g "^ — with 

d = (1.308 ± 0.028) GeV and e = (0.273 + 0.008) GeV" 1 . 
Further the ratio of baryon to strangeness chemical potential is 
^ ~ 0.164 + 0.018 JsZZ. We have checked 



parameterized as — 
that the value of T and p B obtained using the yields extrapo- 
lated to An or from midrapidity measurements, have little im- 
pact on our study. However in order to study the rapidity (y) 
dependence, the p B parameterizations p B - 0.024 + 0.01 ly 2 
and p B = 0.237 + 0.01 ly 2 at yT^ = 200 GeV d and 
■\/s NfJ - 17.3 GeV 11811 . respectively are used in the calcula- 
tions. 



3. Results 

3.1. Kinematic acceptance in r\ and pj 

Figure Q] shows the variation of Xx iXx an d Xx IXx 
as a function of 77 acceptance for y/s NN = 200 GeV and 
■\Js NN = 17.3 GeV. Where x stands for either net-baryon (B) 
(Fig. [2(a) and (b)), net-charge (0 (Fig. [2(c) and (d)), or net- 
strangeness (S) (Fig. [2(e) and (f)). For each beam energy we 
show the effect of considering HRG parameters (p, T) fixed 
to parameterization based on midrapidity data and also a pa- 
rameterization based on the 77 dependent value of (p, T). The 
difference between the two cases are small. For the subsequent 
studies we only present results for different y/s^N using the pa- 
rameterization of the chemical freeze-out parameters based on 
the measurement of particle yields at midrapidity. A clear de- 
pendence of Xx iXx an d Xx IXx on V acceptance is observed 
for net-charge (Fig. [2(c) and (d)) and net-strangeness (Fig. [2 
(e) and (f)). Th &X flxf aadxf/xf values (Fig. [2 (a) and 
(b)) are however found to be independent of 77 acceptance for 
the two beam energies studied. This underscores the need to 
carefully consider 77 acceptance effects when comparing HRG 
model results to experimental data, especially for net-charge 
and net-strangeness fluctuation measures. 

Figure [2] shows the variation of x\ IXx an d Xx IX x as 
a function of V^nn for various pj acceptances. The choice 
of these particular values of pj acceptance ranges are moti- 
vated by existence of the corresponding experimental measure- 
ments J5I [l2j, [l3[] . It is observed that^ v Ixk an d^ r IXx have 
a clear pj acceptance dependence at all beam energies for net- 
charge (Fig.[2](c) and (d)) and net-strangeness (Fig.[2](e) and 
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Figure 1: The r\ acceptance dependence of Xx IXx &®&Xx IXx f° r two dif- 
ferent beam energies. In panel (a) and (b) x = net-baryon B, (c) and (d) x = 
net-charge Q, and in (e) and (f) x = net-strangeness S. Also shown are the 
results with (labeled "dep.") and without (labeled "fixed") the r\ dependence of 
chemical freeze-out parameters yig and T. 



Figure 2: (Color online) The pr acceptance dependence of Xx iXx an d 
Xx IXx f° r different VsnS . Where x stands for either net-baryon (B) (pan- 
els (a) and (b)), net-charge (Q) (panels (c) and (d)), and net-strangeness (S) 
(panels (e) and (f)). 
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Figure 3: (Color online) The variation of x\ iXx an &Xx IXx f° r net-baryon 
(B), net-charge (Q), and net-strangeness (S) as a function of collision energy 
( V- v nn)- The results are shown for different baryons (panels (a) and (b)), elec- 
tric charge states (panels (c) and (d)) and strangeness number (panels (e) and 
(f)) considered in the calculation. 



(f)). However the pj acceptance dependences for net-baryon 
(Fig. 0(a) and (b)) is substantially weaker. Hence the pj ac- 
ceptance study also emphasizes the need to consider the actual 
experimental acceptance for model comparisons in fluctuation 
measures. At the same time both the kinematic acceptance stud- 
ies in r\ and pj shows net-baryon fluctuation measures are least 
affected. 

3.2. Conserved charge states 

Figure [3] shows the variation of Xx iXx an & Xx IXx as a 
function of y«NN for various types of baryons (Fig. |3](a) and 
(b)), values of electric charge states, \Q\ = 1 and \Q\ > 1 (Fig. [3] 
(c) and (d)), and values of strangeness number, |5| = 1 and \S\ 
> 1 (Fig. |3](e) and (f)). For each of the cases the observables 
are compared to the respective values with inclusion of all con- 
served charge states and baryons. We find a strong dependence 

of the x'q /Xn aa ^ Xn 1X0 on wnemer we consider \Q\ = 1 or 
\Q\ > 1, both differing from the case of inclusion of all charge 
states. Same is the situation for net-strangeness. On the other 
hand, successive inclusion of different baryons, starting with 



protons seems to have some small effect on the x B IXb anc ' 
X B iXb values only at the lower beam energies. The absence of 



Figure 4: (Color online) The variation ofXxIXx ar >d Xx iXx f° r net-baryon 
(B), net-charge (Q), and net-strangeness (5) with different beam energies 
( V^'nn) with and without resonance decay daughter particle acceptance effects. 



baryons with \B\ > 1 makes the net-baryon number fluctuations 
more advantageous and less prone to kinematic acceptances as 
compared to the net-charge or net-strangeness fluctuation mea- 
sures. 



3.3. Resonance decay 



In a realistic experimental condition, it is likely to happen 
that not all the daughter particles of a resonance fall within a 
given experimental acceptance. This would lead to a differ- 
ent kind of efficiency factor in determining the net-charge, net- 
strangeness and net-baryon number distributions. For example 
consider a A ++ which decays to proton and pion. Three possi- 
bilities may happen, both the proton and pion could fall within 
the experimental rj and pj acceptances, one of them could fall 
into the acceptance or none of them could fall into the accep- 
tance. These cases will contribute differently to the correspond- 
ing measured net-charge and net-proton distributions. 

We have studied this effect in the following manner. The 
conventional way to incorporate resonance decay is to use the 
definition 111911 . 



(Nd = <jv;> + Y/FxXn)*' 



(5) 



where N* is the primordial yield of the /''' hadron species, Nr 
is the contribution from resonance decay, the summation £r 
runs over all types of resonances and («,)« = £,. bf« R , with 
n,> being the number of particles contributing to the i' b species 
via the decay mode r of the resonance R and b R r is the corre- 
sponding branching ratio. The < «, > can be interpreted as 
a type of efficiency (eg) correction to the contribution coming 
from the resonance decay. Where e« = 1 in the absence of 
any kinematic cuts. The particle yields then can be written as 
(N) = J]i(N)i + J]r £r(Nr). Assuming that the resonance yields 
follow a Poisson distribution, the mean of the distribution can 
be corrected by factor e« and the same correction factor can 
be applied for rest of the cummulants. If N is the number of 
resonances and n is the number that survives the experimen- 
tal acceptance, then e« = bfn/xN, where x is the conserved 
charge - B, Q or S . We have used e« like an efficiency cor- 
rection factor while computing the susceptibilities as given in 
equations |3] and |4] The results of x\ Ixk and^, Ixk with 
and without considering the effect of resonance daughter parti- 
cle acceptances are shown in Fig. [4] In Fig. Ufa) and (b), three 
cases ioT x B /x B an d^ /X B are snown within a realistic ac- 
ceptance of |?7| < 0.5 and 0.4 < p T < 0.8 GeV/c. The results 
for all baryons without any resonance decay (solid red curve), 
results for protons without any resonance decay contribution 
(dotted blue curve) and results for protons with resonance de- 
cay (dashed green curve). Similarly in Fig. |4jc) and (d) shows 
Xn IXn & n d Xn IXn respectively, for all charges without res- 
onance decay (solid red curve), pion, kaon and protons with- 
out resonance decay (dotted blue curve) and pion, kaon and 
protons with resonance decay (dashed green curve). Fig. |4}e) 
and (f), are shown the results for all strangeness without reso- 
nance decay (solid red curve), kaons without resonance decay 
(dotted blue curve) and kaons with resonance decay (dashed 
green curve). As one can see from the figure, the resonance de- 
cays certainly affect all of the ratios. The effects of acceptance 
of resonance decay daughters are large for net-charge and net- 
strangeness compared to net-baryons. The acceptance used in 
this analysis are modest and are close to the present experimen- 
tal scenario. 

Through our work we have emphasized the need for consid- 
ering experimental acceptances of various kinds in model, such 
as HRG, before they are considered to provide the baseline to 
experimental measurements for drawing physics conclusions. 
Hence in the Tables[T][2l[3jand H]we provide values oiXx iXx 
and^ ( v Ixk f or typical ongoing experimental acceptances. The 
Xq IXn &n ^Xn IXn values are provided for two typical accep- 
tances \tj\ < 0.5, 0.2 <p T < 2.0 GeV/c and \tj\ < 0.35, 0.3 < p T 
< 1.0 GeV/c (Tableland Table©. Thexf/xf and^ 4) /*f 
are provided for a typical acceptance of \t]\ < 0.5, 0.2 < pr < 2.0 
GeV/c (Table |3}. Thex B /Xb an d;^ IXb are provided for a 
typical acceptance of \tj\ < 0.5, 0.4 < p T < 0.8 GeV/c (Table©. 



4. Summary 

In summary, using a hadron resonance gas model we have 
studied the effect of limited experimental acceptance on observ- 
ables like n' h order susceptibilities x" > associated with con- 



Table 1 : Ratios of the moments for net-charge within |j/| < 0.5, and 0.2 < pj 

< 2.0 GeV/c 



y/SNN 


Xq IXq 


*M 


5 


0.526 


1.413 


7.7 


0.414 


1.430 


11.5 


0.321 


1.407 


15 


0.265 


1.390 


19.6 


0.215 


1.375 


27 


0.165 


1.361 


39 


0.119 


1.352 


62.4 


0.077 


1.346 


130 


0.038 


1.343 


200 


0.025 


1.342 


2760 


0.002 


1.341 



Table 2: Ratios of the moments for net-charge within |;/| < 0.35, and 0.3 < 

p T < 1.0 GeV/c 



y/^NN 


xf/xf 


*M 


5 


0.432 


1.231 


7.7 


0.332 


1.245 


11.5 


0.256 


1.233 


15 


0.212 


1.222 


19.6 


0.172 


1.213 


27 


0.132 


1.205 


39 


0.095 


1.199 


62.4 


0.062 


1.196 


130 


0.031 


1.194 


200 


0.020 


1.193 


2760 


0.001 


1.193 



Table 3: Ratios of the moments for net-Kaon within |;/| < 0.5, and 0.2 < pr 

< 2.0 GeV/c 



■\[snn 


*M 


xf/xf 


5 


0.726 


1.058 


7.7 


0.560 


1.090 


11.5 


0.428 


1.108 


15 


0.353 


1.115 


19.6 


0.288 


1.120 


27 


0.222 


1.123 


39 


0.162 


1.125 


62.4 


0.107 


1.127 


130 


0.054 


1.128 


200 


0.035 


1.128 


2760 


0.003 


1.128 



Table 4: Ratios of the moments for net-proton within \rj\ < 0.5, and 0.4 < pr 

< 0.8 GeV/c. 



■\[snn 


Xp IXp 


Xp IXp 


5 


0.981 


0.961 


7.7 


0.987 


0.975 


11.5 


0.962 


0.984 


15 


0.920 


0.988 


19.6 


0.851 


0.991 


27 


0.737 


0.993 


39 


0.586 


0.995 


62.4 


0.407 


0.996 


130 


0.210 


0.997 


200 


0.139 


0.997 


2760 


0.010 


0.997 
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served quantities like net-charge (x = Q), net-strangeness (x = 
S) and net-baryon number (x = E). The various order suscep- 
tibilities which can also be calculated in QCD based models 
are related to the moments (<x, S and k) of the corresponding 
measured conserved number distributions. These observables 
have been widely used to understand the freeze-out conditions 
in heavy-ion collisions and various aspects of the phase struc- 
ture of the QCD phase diagram. Our study demonstrates the 
importance of considering experimental acceptances of differ- 
ent kinds before measurements are compared to theoretical cal- 
culations, specifically in the use of HRG model as a baseline 
for such fluctuation based study. We observe finite kinematic 
acceptances in 77 and pj have a strong effect on the^ ' and^J, 
values. These susceptibilities are also very sensitive to the ac- 
cepted electric charge states and strangeness states in the exper- 
iment. In addition we find in this model the;^' and^ n) values 
depends on the experimental acceptance of the decay daugh- 
ters from various resonances produced in high energy heavy- 
ion collisions. Within this model and the kinematic regions 
used in the our study, we find that the dependence on accep- 
tance and resonance decays are stronger for both net-charge and 
net-strangeness compared to that of net-baryons. 
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